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Abstract. The free boundary problem of a two phase flow in a rotating Helc-Shaw cell with 
Coriolis effects is studied. Existence and uniqueness of solutions near spheres is established, and 
the asymptotic stability and instability of the trivial solution is characterized in dependence on 
the fluid densities. 



1. Introduction 

The motion of one or two fluids confined to a narrow gap between two parallel plates is an 
interesting problem with a long history. It is the classical set up of the so-called Hele-Shaw 
problem |13j . It is well-known that instability driven pattern formation such as fingering can occur 
under appropriate assumptions on the viscosity of the fluids |15| and on the surface tension at their 
interface. There is a vast empirical and theoretical literature concerning this classical problem. 
From the purely mathematical point of view the problem has been extensively studied in its original 
formulation as a one [HI El |H1 [JTJ] and as a two-phase QI5] problem and, more recently, also for one 
fluid in the case of rotating plates [7[. The focus is here on the two-phase problem with rotating 
plates and including the effects of Coriolis forces. This setting has recently been considered in the 
physical literature by a variety of authors [TSJ [TTJ [T] . It first appeared in [TB] where the effects of 
rotation were introduced in a ad-hoc fashion into Darcy's law for the one-phase problem 

_ 12?^ 2 _ „ _ _ 

Vp = — -TT-v + quj x + 2guiz x v , 

where z is the axis of rotation, b is the plate spacing, x is the two-dimensional position vector, tu is 
the (constant) angular velocity of the plates, whereas v, p, rj, and g are the fluid's velocity, pressure, 
viscosity, and density, respectively. The last term accounts for Coriolis' force. While most of the 
literature hitherto neglected this force, recent studies performed a model derivation for a one-phase 
Helc-Shaw type model with Coriolis force by means of a standard gap averaging technique starting 
from Navier-Stokes' equations, cf. [TT1 |T5] . The authors of [JS] observe that the effects due to the 
Coriolis term in their equations can be larger than the inertial terms typically neglected in small 
Reynolds number type reductions. They also point out the fact that simpler ad-hoc models used 
earlier, while qualitatively similar, do eventually lead to a different and inaccurate prediction of 
the growth rate of the unstable modes involved in the fingering phenomenon. 

More recently, the effects of a Coriolis term on the fingering patterns in the rotating two- 
phase Hele-Shaw problem have been studied in [TTJ [JJ where the formal linear stability analyses 
of [TBI HH1 for the one-phase problem are extended to the weakly nonlinear case by the use of 
formal expansions and to the fully nonlinear regime by means of numerical simulations. Previous 
numerical and experimental results included |17| 0] , whereas the practical relevance of the problem 
is attested by a number of publications cited in [TJ for that very purpose. 

In this paper local existence of a unique classical solution for nearly circular initial interfaces 
is established for the general rotating two-phase Helc-Shaw problem with Coriolis effects in the 
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formulation proposed in [111 [TJ, which is, in turn, based on a derivation similar to that of |18) for 
the single phase case. It leads to the following generalized Darcy's law 

for the velocities of the inner and outer fluids (j = i, o), for constants ay and (3j defined below and 
for the pressure Pj related to the hydrostatic pressure pj via 

The use of the unusual pressure Pj yields some useful simplifications as will soon become apparent. 
Rigorous stability results are also obtained on near circularity assumptions for the initial data by 
resorting to a general principle of linearized stability. While the circular steady state is exponen- 
tially asymptotically stablqj when the outer fluid is denser, it becomes unstable if this relation is 
reversed. The main techniques used are a transformation of the original free boundary problem to 
a fixed domain one, a decoupling of the system and a reduction to a single nonlinear and nonlocal 
evolution for the interface separating the two fluids, and parabolic optimal regularity results in 
little Holder spaces. 

Related previous mathematical results include the existence, stability, and bifurcation analysis 
of [7] in the rotating one-phase problem and the global existence and stability of smooth solutions 
obtained by |19j for nearly circular initial interfaces in the non-rotating two-phase case. 

2. Governing Equations and Main Results 

Wc first give a short justification for the mathematical formulation of the physical problem of a 
rotating Hele-Shaw cell including Coriolis effects. Further details regarding the modeling can be 
found in [TJ [TTJ . We then end this section with the main results on existence of solutions to the 
corresponding governing equations and stability properties of the trivial solution. 

2.1. Governing Equations. Consider a circular Hele-Shaw cell of radius R > 2 and very small 
gap width b > 0, rotating clockwise around the z-axis with a constant angular velocity u > 0. The 
cell is assumed to contain two immiscible, incompressible, viscous fluids with densities Qj > and 
viscosities r\j > 0, where j =i labels the inner and j —o the outer fluid, respectively. The surface 
tension between the two fluids is given by a > 0. The rotating coordinate system is defined in such 
a way that its origin is located at the cell center and that rotation is perpendicular to the vector 
z = (0, 0, 1). Let f2-?'(i) denote the region of space occupied by fluid j £ {i,o} at time t, and let 
r(i) denote the sharp interface separating the two fluids. The unit normal vector ^r(t) to T(t) is 
assumed to point from fi'(i) to fi°(t). 

The basic hydrodynamic equation of the system is a generalized Darcy 's law relating the pressure 
fields Pj = Pj (t) to the two-dimensional gap-averaged flow velocities Vj = Vj (t) through 

VPj = -oijVj + Pj(zx Vj) in Cl j (t) . (2.1) 

The numbers 

a 3 ■ & ^ ' Pi — 1)2 P 3 

depend linearly on the Coriolis force terms Ej > and Fj > 0, respectively. The latter, in turn, 
depend on the rotational Reynolds number Rej = QjUjb 2 /12rjj (see [TJ). If Coriolis effects are 
neglected, then Ej = 1 and Fj = 0, so (|2.ip reduces to the usual Darcy's law. Throughout this 
paper we shall assume that ay > and j3j > 0. Incomprcssibility of the fluids is expressed by 

divu i = in Cl j (t), (2.2) 



-^Our analytical approach reveals a precise description of the exponential decay in terms of the physical param- 
eters, in particular of the Coriolis terms, see Section [4] 
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while the interface dynamics is governed by the normal stress balance 

Pi-P = <7rt r(t) + (7 -7i)|z| 2 on r(t), (2.3) 

and the kinematic boundary condition 

V = n ■ vv = v ■ vr on T{t) . (2.4) 

Here, jj := QjLo 2 /2 whereas Kr(t) denotes the curvature of T(t) which is taken to be positive if 
Q l (t) is convex. The normal velocity of the interface is given by V — V(t). On the outer boundary 
of the cell a no-slip condition 

v Q ■ x = on [|x| = R] (2.5) 

is imposed. 

Observe that the inner and outer fluids' masses are conserved since, by Reynold's transport 
theorem, one has that 

d r f _ . r ^ x , 

ax = — / v Q ■ Vr dcr + / v Q ■ — da 



Jn°(t) Jr(t) •AkN-R] ^ 

and 

— / dx = I vi-vrda, 
at Jsi l (t) Jr(t) 

which are both zero in view of Gauss' Theorem and (|2.2p . Since rot(z x Vj) = (0, 0, diviJj) we also 
note that the flows are irrotational in the bulk, i.e. notvj = in f2 J '(t), due to (|2.1j) and (12. 2ft . 

2.2. The System in Terms of Pj. The governing equations in terms of Pj are obtained by 
taking the divergence on both sides of (|2.1[) which yields 

AP,=0 in n j (t), 

making use of (|2.2I) and observing that rotiJj = implies div(z x Vj) — 0. As for the kinematic 
boundary condition (|2.4p . we solve (|2.1[) for Vj to get 

1 



|0 3 f 



with complex numbers 0j := q_, + Then, defining the tangent vector rp := — (z x to T, 
we obtain 

^ ' ^ r = |0l2 ( ~ a 3 d vr P 3 - Pj d Tr P j) ■ 

Whence (l2~4l) and (1231) become 



-V = (aid^Pi + ASrpPi) = -J— (ttoft^Po + Mr r Po) On T(i) , 

|6ir |© r 

and, 

a<AP + &AP = on [\x\=R], 
respectively, where v denotes the unit outer normal vector and r the corresponding tangential 
vector on the cell boundary [|x| = R\. Therefore, we arrive at the following free boundary problem 
for the pressures Pj 

APj=0 mfl j (t), J=i,o, (2.6) 

P i -P = an m + ( 7o - 7i )|x| 2 on T(t) , (2.7) 

-V = — ^—(aid^Pi + A9 rr Pi) = — ^(Oo&rPa + MrrPo) on T(t) , (2.8) 

|©iT |6 | 

a<AP + p o drPo = on = R] , (2.9) 
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for t > , complemented with an initial surface T(0) = T° . Since only derivatives and the difference 
of the two pressures enter the system, uniqueness can only be expected up to additive constants 
for them. As we shall later see, this is one of the obstacles that need to be overcome from a 
mathematical viewpoint. 

2.3. Main Theorems. To give a precise formulation of our mathematical results on (|2.6[) - (l2.9p . 

we parametrize the boundary T(t) over the unit sphere S 1 := {x G R 2 ; |x| = 1}. To this end 
we introduce, for s > 0, so-called little Holder spaces h s (U) over an open subset U of R n as the 
closure of BUC^iU) in BUC S (U). Here BUC a (U) consists of all functions / : U -> M with 
bounded and uniformly continuous derivatives up to order [s] and with uniformly (s — [s])-H61dcr 
continuous derivatives of order [s]. If M is a (sufficiently smooth) submanifold of R™, we define 
h s (M) by means of an atlas for M in the canonical way. In the following we shall identify a 
function h : S 1 -> K with the function h : [0, 2tt] -> M given by h(6) := h(e w ) for < 6 < 2ir. The 
bar will often be dropped as no confusion seems likely. 
Fix a 6 (0, 1/4) and 8 G (0, 1) and set 

V^ipGh^iS 1 ); ||p|| 00 <o} l 

and, for p G V, 

nj,:= {xGK 2 \{0}; |*| <l + p(^)}u{0}, 
with ft := B(0, R) where, we recall, i? > 2. Then 

T p := L G M 2 ; |z| = 1 + p | = { [1 + p(y)]y; y G S 1 } 

separates the domains ftp and ft°. Since T p can be described as the zero level set of 

N„(x) := \x\ — 1 — p ( j~t) j 3/4 < |x| < 5/4, 



with N p < in ft' n [3/4 < |x| < 5/4], the unit normal ^ p (x) at x G T p pointing from ft p to ft° is 
given by 



In the following, we let r p = — (z x i/ p ) denote the corresponding tangential vector. Next, suppose 
that p G C([0,T],V) nC 1 ([0,T],/i 1 + f (S 1 )) for some T > and set 



JVpftar) :=JV p(t) (ar) = |z| - 1 - p ^, ^ , te[0,T], 3/4 < |x| 



< 5/4 



Observe that the normal velocity V p of the moving boundary T p equals — dtN p /\VN p \. Conse- 
quently, if p G C([0,T], V) n ^([OjT],/! 1 ^^ 1 )) describes the evolution of the boundary, then 
(1211)- (Ell]) read as 

AP J= i = V>, (2.H) 

P» - Po = <r«r pW +(7o-7i)[l + p(t^)] 2 on T p(t) , (2.12) 

M f > o) = jeqi h^Pj + x VP;)] • V7V p(t) on T p(t) , j = i,o , (2.13) 

a d v P + MrPo = on ,90 = PS 1 , (2.14) 
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for t > and with r p / ) = T° . We call a triple (p, Pi,P ) a (classical Holder) solution to (|2.11|l - 
dHIl| provided 

p £ C([0,T],V) nC 1 ([0,r] J /i 1+ '(S 1 )) 

for some T > and 

i , i(t,-)G/» a+ '(nJ w ), ie[0,T], j = i,o, 

satisfy (|2.11l) - (|2.14[) pointwise. The main existence result is stated in the following theorem, a 
proof of which is given in Section [3J 

Theorem 2.1. Let S £ (0, 1). There is an open zero neighborhood O in /i 4+5 (§ 1 ) such that for each 
initial geometry pa £ O, there exist a time T := T(po) > and a classical solution (p, Pi,P a ) to 
(I2.11[) - ([2.14j) with p(0) = po- This solution is unique except for additive constants in the pressures 
Pi and P . 

We also investigate stability properties of the trivial solution corresponding to the geometry 
p = 0, that is, to the unit circle S , and constant pressures P\ = c + a + j Q — 7; and P = c with 
arbitrary c £ R. The stability properties of the trivial solution are determined by the relative 
density of the fluids. 

Theorem 2.2. If g Q > q\ then the trivial solution (p,Pi,P ) = (0, c + a + 70 — 7i,c) is locally 
asymptotically stable. However, if > g , then the trivial solution (p, Pi, P a ) = (0, c+ct+7 — ji, c) 
is unstable. 

We refer to Section [4] for more precise statements of the stability results and their proof. 

3. Proof of Theorem 12.11 

The proof of Theorem 12.11 is best carried out in a coordinate system in which the moving 
interface between the liquids is fixed. We therefore begin with the transformation to a fixed 
domain formulation. 

3.1. An Equivalent Problem on Fixed Domains. We transform the free boundary problem 
to fixed domains using the standard Hanzawa-transform. Let ip : R — » [0, 1] be a smooth function 
with tp(r) = 1 for \r\ < a and f(r) = for r > 3a and ||</s'||oo < \ja. We introduce a family of 
C 4+(5 -diffcomorphisms 

4> p £ Diff 4+5 (R 2 ,R 2 ) nDiff^sSnj,) nDiS 4+5 (B°,n°) , peV, 

where 

B i —3(0,1), B°:=Cl\B 1 , 

by setting 

0p(x):= /(N + [N-iHr))r ifo<M<2, 

[x else. 
Note that 4> p maps S 1 onto T p . Given p £ V, let 

4>* p : BUC(W p ) BUC(B* ), p^ P o$ p 
denote the push-forward operator and 

4>p : BUC(Bi) ->• BUC(QP p ) , q^qo^ 1 

the pull-back operator induced by <f> p . Given p £ V, the transformed differential and boundary 
operators acting on Qj := <fi p Pj j = i, o are given by 

Mp)Qj :=r P (H€Qj)) 
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and by 

B MQi [«i(^v(0SQ 3 -)) + fr{zx (#v(#Q 3 -)))] • (^;viv p ) . 

Defining 

K{p) :=a^ Krp + (7o-7i)|l+p| 2 

for p G V and 

B := a Q d v + j3 d T on 90 = RS 1 , 
the free boundary problem (|2.1ip - (|2.14[) is transformed to the following problem on fixed domains: 



Aj(p)Qj = inB\ j = i,o, (3.1) 

Q l -Q = K{p) on S 1 , (3.2) 

d t p = Bj(p)Qj on S 1 , j = i,o , (3.3) 

B Q = on RS 1 . (3.4) 



We call a triple (p,Qi,Qo) a (classical Holder) solution to (|3.1[) - (|3.4[) provided 

p e C([0, T], V) n C 1 ([0, T], ^(S 1 )) 

for some T > and 

Q. CV) e /j 2+5 (B 3 ') , te[0,T], i = i,o, 

satisfy (pTTjl - tpO)) pointwise. With this definition, problems ([2~TT|) - (|2~T4]) and ([3TT ]) -(|3~4" |) are equiv- 
alent as follows from the next proposition. 

Proposition 3.1. Let Qj = <p*Pj and Pj = $Qj. Then {p,P h P Q ) is a solution to (|2.11|) - (|2.14p 
if and only if (p, Qi, Q Q ) is a solution to (|3.1[) - (|3.4[) . 

Proof. Due to the above definitions of the differential operators, we merely need to ascertain that 
no regularity loss occurs in the process, i.e., that Pj G h 2+s (tti p ) implies Q 3 = <j>*Pj e h 2+s (B^) 
for p G V and vice versa. This, however, easily follows from the fact that <fi* p Pj G BUC A {B : ') 
when Pj G BUC°°{tt j p ) and p G V C /i 4+5 (§ 1 ) and the observation that h 2+s {B^) coincides with 
the closure of BUC 4 (B^ ) in BUC 2+S (B^). □ 

Next we collect some regularity properties of the operators defined above. 

Lemma 3.2. (i) The operator Aj(p) G C(h 2+s (B :l ),h s (B :l )\ is uniformly elliptic and analytic in 
peV for j = i,o. 

(ii) The operator Bj{p) G C(h 2+ (B^),h (S 1 )) is analytic in p G V /or j = i,o. 

Proof, (i) Given p G V, uniform ellipticity of -4j(p) is a consequence of its symbol being 

= j(0 P - 1 )ei 2 + • e) , e g k 2 , 

where the Jacobian J((j>~ 1 ) of 0" 1 satisfies |J(</>p 1 )C| 2 > e l£| 2 f° r some e > 0. Analyticity in p 
follows from the analyticity of <f> p and of J(4> p 1 ) = (j)*J(4> p ). For details we refer to [H § 3.2]. 
(ii) Note that 

VP (r) =P '(r)^(- 2 )' *=(^)e* a W°>. P^W. (3.5) 
Hence, for p 6 V and 4> p = (4> p , </> 2 ), 

u , ^(y) t ( <t> P {y) \ i {-<i>e(y)\ 
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showing that [p \-¥ </>*ViVp] , V — > /i 3+5 (§ 1 ) is analytic. The definition of Bj(p) entails its analytic 
dependence on p since 



□ 



The curvature operator K, can easily be computed explicitly 
Lemma 3.3. The operator K. : V — > h 2+s (S 1 ) is analytic and given by 

MW = ff " + (7o -7i)(l +p) , /)£V ; 

[(l + p) 2 +p 2 ] 

and 

5/C(0)[/i] = a(-h -h) + 2( 7o - 7i )/i, h e ^(S 1 ) . 
Proof. If p <E V, a parametrization of T p is given by 

(1 + /J(#))e* e =: x{6) + iy(0) , G [0, 2tt] . 
Hence the claim is a consequence of the curvature formula 

xy — xy 



. 2 n3/2 

y 2 ) 



and the definition of K,(p). □ 

3.2. Local Well-Posedness. In order to establish local existence and uniqueness of solutions, we 
split p.lj) - (|3.4[) in three subproblems: two elliptic problems for the transformed pressures and an 
evolution equation for the geometry. More precisely, given p G V and a function h defined on S , 
we shall first look for solutions Qi = S(p, h) to the following elliptic problem on B 1 with Neumann 
type boundary condition 

Ai{p)Qi =0 in B l , (3.6) 

Qida= [ hda, (3.8) 

where |r p | stands for the length of the closed curve T p . In a second step, we shall study the elliptic 
problem on the annulus B° with Dirichlet-Neumann boundary conditions 

Mp)Qo = inB°, (3.9) 
Qo = g on S 1 , (3.10) 
B o Q o = onfiS 1 , (3.11) 

with solution Q = T(p,g) depending on p S V and g on S 1 . Finally we use the solution operators 
S and T to derive the evolution equation 

d t p - B (p)T(p, S(p, d t p) - K{p)) . (3.12) 

for p. Note that dtp appears on both sides of (|3.12[) . As (|3.6[) - (|3.7[) is a Neumann- type problem, an 
additional integral term is introduced on the right hand of side (|3.7[) which makes the mean zero 
and thus ensures solvability of the problem. The third equation (13.81) guarantees uniqueness. Let 
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us point out that, provided we can solve (|3.12[) for p = p(t), the particular form of the boundary 
operator B (p) and the properties of T(p,g) yield (see (|3.44p below) 

Jt p |vj\g 

Thus the integral term in (|3.7I) vanishes for h = dtp and we may take 

Qi(t) = S(p(t), d tP (t)) , Q„(t) = T(p(t), Qi(t) - K(p(t))) 

to obtain a solution for the original equations (|3.1I) - (|3.4[) . We shall be more specific at the end of 
this section. 

In the following two propositions we study the solvability of the elliptic problems. 
Proposition 3.4. Given p £ V and h £ h 1+d (S 1 ), (|3.6I) - (|3.8I) possesses a unique solution 

Qi = S(p, h) £ h 2+s {B l ) , 

and the map 

\pH>S(p,-)] : V-^^+'tS 1 ),/* 2 ^ 1 )) 

is analytic. 

To prove this statement on the fixed domain L? 1 we shall first prove an auxiliary result on the 
domain f2* : 

Lemma 3.5. Let p £ V and Ze£ i> = if and t = 7r denote the corresponding outer unit normal 
and unit tangential vectors to T p , respectively. Set 

and define 

A : h 2+5 (Q p ) -> h 5 (fl p ) x h 1+5 (T p ) , ui->(-Au, d^u) . 
T/ien kcr(A) = R • 1 and 

im(A) = |(/, 3 )£/ l 5 (^)x^ 1 + 5 (r p );^ fdx=^J^ gdaj . 

Proof. Note that is nowhere tangential and that 1 is an eigenfunction of A. Thus, kcr(A) = K • 1 
follows from Thm.12.1] since (— A,^) is regular elliptic. To determine the range of A, fix 
p > n/(l — S) and suppose first that Au = (/,<?), that is, 

- Ait = / in Cl p , d^u = g on T p . (3.13) 

Since d T w = (z X Vw) • 1/ and div(z x Vw) = 0, we have by Gauss' Theorem 

/ d T wdcr= I div(z x Vw)dx = 0, w £ W 2 (nj,) . (3.14) 

Therefore, 

/" , , |e ; | 2 /■ / ai _ A _ A , |6i| 2 /■ 

/ /dx= / -ttvti^u- -— <9 T n do-= / gda . 

Ja p «i Jr p V l il l©il 2 / "i Jr p 

For the reverse inclusion we use a Frcdholm argument. By [21 Lem.5.1] there is a coretraction 

M££(w p 1 -^(r p ),w; 2 (^)) 

for the boundary operator such that d p Mg = g for g £ W^ 1 1 ^ p (r p ). Then 

i?:=Ai\/£/:(^ 1 - 1 / P(rp);jLp(fi i ) ) i 
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and finding a solution u e Wp(Q p ) to problem p,13[) for a given (f,g) G £ p (fip) x Wp -1 ' p (r p ) is 
equivalent to finding a solution v £ Wp(Cl p ) to 

-Av = f + Rg in , <V> = on T p , (3.15) 

and setting u := u + Mg. Let 

% 2 p (r> p ) := {v G W^ P ) i = on T,} 

and consider the closed linear operator T on L p (fl p ) given by Tv := —Av for v € Wpu(^p)' Then, 
as above, ker(T) = K • 1, and T has compact resolvent. Hence T is a Frcdholm operator on L p (£l p ) 
with index zero. Since obviously W p (fl 1 ) is dense in L p (il p ), we also obtain for its dual T' that 
ker(T') = M • /' for some /' € Zy (O p ), where p' denotes the dual exponent of p. Moreover, Tv = h 
with h G L P (Q„) is solvable for w G W p JCl p ) if and only if (f',h) = 0. However, since (|3.14p 
ensures 

/" A vdx = - l -^[ d p vda = 0, veW^(n [ p ), 
Jn p oti j Vp 

we may take /' = 1, that is, ker(T') = K • 1, and Tv — f + Rg is solvable for v G W pp (f2 p ) with 

(f,g)GL p (n p )xW^(T p ) 

if and only if 

= / (/ + Rg) dx = ( fdx- / 5 p .g da , (3.16) 

Jcfi„ Jcfi p on J Fp 

the last equality being again due to (|3.14|) . Finally, if u is the solution to ()3.13j) for a given 

(/,«?) g h s (n p ) x /. 1+5 (r p )cL p (fi P ) x w^(r p ) 

satisfying (|3.16p . then u G Wp(fl p ) <^-> C 1+(5 (f2 p ) by the choice of p. Using Schauder regularity 
theory [T3J Thm.6.15] for 

— Ait = / in Q 1 , u + <9 p u = g + u on T p , 

then easily gives u = u G /i 2+<5 (f2 p ). This proves the claimed characterization of im(A). □ 

Proof of Proposition \3.4\ For p G V fixed, set 



w 'Tfl,ra'" " e " <s> ' 

and define an operator A i (p) : h 2 + 8 {B') ^(B 1 ) x ^ 1+,5 (S 1 ) by 

Ai(p)u := (fl(/))A(p)«,(!(p)fii(/))ti), u G /i 2+<5 (_B i ) , 

where -D(jo) := |det 9^0^,1 and d(p) := 1 + p. Then, since D(p) > and d(p) > 0, we readily obtain 
from Lemma l3.5l that ker(Aj(p)j =R-1. Moreover, (f,g) G im(A;(p)) if and only if (/, g) G im(A), 
where 

By virtue of p.5[) we have that 



|V7V p |= VP+(1 + P) 2 , on sl (3 . 17) 
1 + p 



and we infer from Lemma 13.51 that 



im(A 5 (p)) = {(/,«?) G h S {B { ) x ^(S 1 ) ; j fdx = f 9^] =■ Y 
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which is independent of p. Consequently, the map 

Ai(p) : ^(B'J/m — > Y 

is an isomorphism. Furthermore, v = Af(p)(f,g) := A- i (p)~ 1 (f,g) is the unique solution in 
h^iB' 1 )/^ to 

D(p)Ai(p)v = f onB 1 , d(p)B i (p)v=g on S 1 
for each p E V and (/, g) 6 F. As and o!(p) are analytic in p, we deduce that 

Af:V^£(Y,h 2+s (B 1 )/^) 
is analytic. Now, we easily check with the help of (|3.17|) that 

d(p)(l - P(p))hda = 0, 

whence (0,d(p)(l - P{p))h) E Y for h E /i 1+5 (S x ). Thus u := AA(/9)(0, d(/t>)(l - P(fi))h) is the 
unique solution in h 2+s (B 1 ) / R to 

A(/»)u = onff, Bi(p)u= h- P{p)h onS 1 . 

By (|3~T7f we have 



C 1 + p) Jo Vp+( 1 + p) d0 Jo 

hence P(p) depends analytically on p and so does u = u(p). Given p E V and h E /i 1+<5 (§ 1 ), 
problem (|3.6[) - (|3,8|l thus admits a unique solution Qi = S*(p, /i) and 

is analytic. This proves Proposition [3^4] □ 



Proposition 3.6. Given p E V and g E ft. 2+5 (§ 1 ), £/iere is a unique solution Q Q = T(p,g) E 
h 2+s (B°) to (I331)- (I3~TU1) . Moreover, tfte map [p h-> T(p, •)] : V ^(^(S 1 ), h 2+s (B°)) is' ana- 
lytic. 

Proof. Let p e V and define an operator A a (p) : h 2+s (B°) -S- x /^(S 1 ) x /i 1+5 (i?S 1 ) by 

A (p)u := (,4o0>)u, u| s i , B u), u E h 2+s (B°) . 

Then A Q (p) is invertible. Indeed, to check its injectivity, let A a (p)u = 0. Then, for u := <j)*u, 

Au^O in fl° , 
u = on , 
a d v u + (3 d T u = on i?§ 1 , 

and Green's formula yields 

uAudx^- [ |Vw| 2 da:+ f ud u uda = - [ |Vu| 2 dx - — [ \d T u 2 da 

| V S| 2 dx - A. y div f dx = - [ |Vu| 2 dx 
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This entails u = 0, i.e. u = 0. To determine the range of A Q (p) consider again the transformed 
problem 

Au = f in Q°, (3.18) 
u = g on r„ , (3.19) 
a a d v u + (3 Q d T u = ft on i?§ 1 . (3.20) 
By means of a coretraction as in the proof of Lcmma l3.5l we may assume g — and ft — 0. However, 

by 0, 

A - A : {<u G W 2 2 (f2°) ; « = on T p , (a Q d v + Mr)v = on RS 1 } — > L 2 (fi°) 

is invertible for A > sufficiently large, and the argument above to prove injcctivity implies 
that we actually may take A = 0. Thus, (|3.18[) - ()3,20[) is uniquely solvable for each / G £2(^0), 
g G iJ 3 / 2 (r p ), and ft G H 1 / 2 (RE> 1 ). Schauder regularity theory ensures w 6 ft 2+,5 (fi°) provided 
/ G h s (ft°), g G ft 2+,5 (r p ), and ft e /^(.RS 1 ). Consequently, A D (p) is surjcctivc and thus 
invertible, and the claim on the analyticity of T(p, •) = A ( ( o) _1 (0, ■, 0) follows from Lemma l3"T2"1 □ 

Next, we state a multiplier result that we shall use later. The proof is a straightforward modi- 
fication of the case r = s in [3] . 

Proposition 3.7. Let r, s G (0, 00) \ N and (M„)„ 6 z be a sequence in C satisfying 

(i) su P«ez\{o} \n\ r ~ s \M n \ < 00, 

(ii) sup„ eZV{0} |?i| r - s+1 |A/„ + i - M„| < 00 , 

(hi) sup„ eZ y {0} |?i| r - s+2 |A/„ + 2 - 2M n+1 +M n \ < 00 . 
Then the mapping 

h n e me .— ► ~ h nMne m0 

belongs to ^(C^S 1 ), C^S 1 )) . 

We now focus on the evolution equation (|3.12[) . which we may rewrite, using linearity of the 
solution operators 5 and T, as 

(1 - K{ P ))d t p = B ( P )T{p, -K(p)) (3.21) 

with 

TZ{p)z:=B (p)T(p,S(p,z)), z G ft 1+5 (S 1 ) , P eV. 
To solve p.21|) for dtp, we use the following 

Lemma 3.8. There is an open zero neighborhood VV C V in h 4+s (E> 1 ) such that 1 — IZ(p) is an 
isomorphism on h 1+s (S 1 ) for each p G W. 



Proof. By the smooth dependence on p stated in Lemma [331 Proposition [33] and Proposition [ 
it suffices to prove that 1 - 11(0) e £(ft 1+5 (S 1 )) is invertible. To do this we derive its Fourier 
expansion in polar coordinates (r, 6). To compute Q; = 5(0, ft) we first note that, for p = 0, 
problem (|3.6I) - (|3.8[) becomes 



\d r {rd r ) + i$ ) Qi = in [r < 1] 



-1 , „ „ „ , „ , 1 



2zr 



^.(a i d r -p i d d )Q i = h-— h(9)d6 on [r = 1] 

l"i| 27r J0 



/ Qi(l,6)d6= / ft(0)d0, 
Jo Jo 
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which, for a given h G h 1+s (S 1 ) with expansion 

h{9) = ^ h„e ine , (3.22) 



has the unique solution 



Q i (r,e) = h + V h n |6i ' 2 , r ^e me . (3.23) 

«eA {0} -H«i + ««A 



Similarly, for /? = 0, problem (|3.9|l - (|3.10|) reads as 
'!„,„, 1 



^ r (rd r ) + ^a|J Q = in [1< r < R] , 

Qo = 9 on [r = 1] , 

a <9 r - Qo = on [r = J?] , 

and for g e ?i 2+<5 (S 1 ) with expansion 

ti6Z 

the unique solution Q Q = T(0,g) is given by 

Oo(r,ff) = A,+ £ 4e^ r " + eTw4 m9 (3 ' 24) 

Therefore, given h E h 1+s (E> 1 ) with expansion (|3.22[) , we have 

n(0)h = (a Q d r - Me) T(Q, 5(0, h)) = £ ^ l ^ mB , 
|Ho1 nez\{0} 

where 

fl — R 2n ) 10- 1 2 

, , { . R ' fn '' R2nh , , «eZ\{0}. (3.25) 

(sign(n)ai - lAj (©o + i? 2n ©o) 

We next use Proposition 13.71 with M n := (1— In)" 1 to check that 1 — 71(0) is invertiblc on C 1+5 (S 1 ). 
Note that l B ^l for each n G Z \ {0} since Oj > 0, j = i, o. Also, 

lim l n = -T^-ff ^ 1 , Um l n = — jT^f 7^ 1 (3-26) 

n->oc |"o| n->-oo |"o| 

so that 

sup |M n | < 00 . (3.27) 

nSZ\{0} 

Next, 

where, for n > 1, 

Q(i? 2 "+ 2 ) OQR 2 ") 

n+1 " " e 2 i? 4 «+ 2 + o(i? 2 ™+ 2 ) ' n - ' n+1 n ~ e 2 + o(i? 2 ™) ' n ~ ' 

Since nR~ 2n —> as n — > 00 and ni? 2 ™ — > as n — > —00, it follows from (|3.27l) that 

sup |n| |M n+ i - M n \ < 00 . 
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Finally, as above, we obtain from (|3.26[) and (|3.27p that 

ln+2 — 2Z Jl+ i + l n + l n (l n +l — ln+2) + ln+2 (ln+1 ~ In) 



n 2 \M n+2 - 2M n+1 + M n \ = n 2 



< (sup|M n |) 3 [(1 + \l n \)n 2 \l n+2 - l n +i\ + (1 + \l n +2\)n 2 \l n+1 - l n \] 

< csup (n 2 \l n +i — ln\) < oo . 

n 

Consequently, Lemma [37F] implies (l - ^(0)) _1 G £(C 1+,5 (§ 1 )). But then (l - ft(O))" 1 is also a 
bounded operator on 

H'(S 1 )=lh€L 2 {S 1 ); \\h\\ H s :=J2(l + n 2 ) S K\ 2 <°4 

I nSZ J 

for each s > due to (|3.27[) . Recalling that -ff^S 1 ) is densely embedded in h 1+s (E> 1 ) provided 
s > 5/2, we deduce that (l - ^(0)) _1 6 C(h 1+5 ($ 1 )) . □ 

According to Lemma 13.81 and (|3.21[) , we are led to look for solutions 

p E C([0,T],W)nC 1 ([0,7 , ],/i 1+< (S 1 )) 

to the fully nonlinear equation 

d t p=(l-K{p))- X B {p)T(p,-K{p))=:F(p), t G (0,T] . (3.28) 

The following proposition is instrumental in the computation of the linearization in zero of this 
evolution equation. 

Proposition 3.9. F G C°° (W , h 1+s (S 1 )) and, for h G ^^(S 1 ), 

dF(0)[h] = (l-^(0)) _1 S o (0)T(0, (2(ti - 7o) + &)h + ah) . 
In particular, if h G /i 4+<5 (§ 1 ) with h{6) = J^nez h n e mf> , then 

dF(0)[h}(9)= ^ne me 

nGZ\{0} 



where, for n G Z \ {0}, 
with 

i? 2n + 1 



A n +isign(n)B 
In ■= A2 +B2 M(n) (3.29) 



An := si 



ig n ( n ) R2n _ 1 a ° + Qi ' B -=Po-Pi, p{n) := |n|(cr-2(7 -7i)-crn 2 ) 



Proof. Smoothness of the map -F follows from Lemma 13.31 Proposition 13.41 Proposition 13.61 and 
Lemma Let /i G /i 4+5 (§ 1 ). Then 

dF(p)[h] = 9(1 - 7e(-)) _1 (0)[/i] Bo(0) T(0 , -/C(0)) + (l - ^(O)) -1 9S o (0)[/i] T(0 , -/C(0)) 

+ (l-ft(0))~ 1 £ o (0)S p T(0, -/C(0))H + (l-^(0)) _1 6 o (0)T(0, -dK(0)[h]) . 

But — /C(0) = — cr — 7o + 7i =: c G M. and T(0,c) = c by uniqueness, so dB o (0)[h] = as this is 
the derivative of (p h->- B (p)c = 0) and similarly d p T(0 , — /C(0)) [/i] = as this is the derivative of 
(p i— > T(p, c) = c). The formula for dF(0)[h] now follows from Lemma l3.3l To compute its Fourier 



14 



JOACHIM ESCHER, PATRICK GUIDOTTI, AND CHRISTOPH WALKER 



expansion, consider h e /i 4+<5 (§ 1 ) with h(9) = E„ez^™ e *" 6 '• Invoking and recalling that 

£> o (0) = — (a d r — (3 de)/\ < d \ 2 on S 1 , we obtain 

dF(O)[h}(0) = (1 - ft(0)) _1 B o (0)T(0 , (2( 7i - 7o ) + a)h + ah) (6) 

_ , n(2( ri ~ 7o ) + a-an 2 )(R 2 "-l) 1 mfj 

J^m e + i?2«e i-i n e 

nGZ\{0} 

with Z„ given by ()3.25j) . Elementary calculations now lead to the assertion. □ 

Observe that A n = A- n for n E Z \ {0} and that 

^4n \ cto + oti as |n| — > oo . (3.30) 

The next proposition is a consequence of the previous lemma and fundamental for our well- 
poscdncss result. 

Proposition 3.10. -&F(0) 6 W(ft 4+5 (S X ), ^ 1+5 (S 1 )) , that is, dF(0) E £(/i 4+ *(S 1 ), /i x+ *(S 1 )) zs 
</ie generator of an analytic semigroup on /i x+<5 (S x ). 

Proof. Based on PropositionEIZl we prove in a first step that (A-aF(O))" 1 £ £(/i x+<5 (S x ), ft 4+<5 (S x )) , 
where, according to Proposition 13. 91 (with qo :— 0), we have that 

(A - 8F(0)) ~ l h = J2 K t - ~ — e in8 , Re A > A* , 



for fc(0) = Y. n& h n e m0 and, by flUO]), 



A* := 1 H > Req n , n E L . 

a a + a; 

Let AeC with Re A > A* be fixed and set := (A — g n ) _1 . Then 



sup n 



3 I r>A I 



(3.31) 



since 



a Q + on ± z_B 

±- 0~ as n > ±00 



(Q + ai) 2 + B 2 
m 2i? 2 »(i?-l) 



Next observe that, for each meN, 

W m l^-^l = H m (i^-i)(i^i-i) 

Thus, letting z„ := q n / fi(n) we derive from (|3.30l) 

A n+1 A, 



as n — > ±oo 



|nr|z " +1 Znl ~ (A 2 +1 +B 2 )L4 2 +B 2 ) 



so, for each m E N, 



Thus, (|3~3T1) gives 



n| m K +1 - A ;l | + \n\ m \A n+1 - A n \ B 2 , 
n\ m \z n+ i — z n \ — >0 as \n\ — > oo . (3.32) 



n | 4 |Qn+l Qn\ — W^Qn+l] \ n3 Qn\ 



Qn+l — Qn 



H(n) 



, \ z n+\\ \ji{n+\) -/i(n)| 
< c - h c z„ + i - z n \ \n\ 

Taking into account that (z n ) is bounded and that |^(n)/n 3 | — >■ a as \n\ — > oo, we conclude 

|4 l^iA 



sup |n| 4 |Q£ +1 -Q£| <oo 

nGZ\{0} 



(3.33) 
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IS 



sup n 2 \q n+1 - q n \ < oo 



In particular, we have shown 



Observe then that 

(i(n) [n(n + 1) - fi(n + 2)] + fi(n + 2) [fx(n + 1) - n(n)] = 0(n 4 ) , 
hence, due to (|3.32[) , 



(3.34) 



1 



9n (Qn+l — Qn+2) + <7n+2 (Qn+l — 9n) 

fi(n) [n(n + 1) - m(" + 2 )] + K n + 2 ) + 1) - /■»(")] 



< 



/i(n) M(n+1) 2 
z« — 5 5 (z n +i ~ z n+2 )n 



fi(n + 2) n(n+l) , r 
z n +i 5 3 \z n +2 — z n ) n 



H(n) n(n + 2) 2 
— s s [Zn+i - z n+2 )n 



(3.35) 



< c . 

Writing 

|^ +2 -2Q^ +1 +^| 



— | n3 Qn+2| I^^Qri+ll | n3 Qn' 

we deduce from (|3.3ip and (|3.34p that 

SUp |n| 5 \Qn+2 ~ 2Qn+l + Qn\ < °° 



A(g n +2 — 2g n+ i + 

(ln+1 — 9n+2) + <?n+2(?n+l — ?n) 



(3.36) 



Consequently, (A - <9F(0)) 1 G £(C 1+J (S 1 ),C* +tf (S 1 )) for Re A > A, by Lemma O and (jXBTD . 
(|3.33[) . and (|3.36[> . Since f|3.31[) also ensures 

(A-9F(0)) _1 g £(F S (S X )) , s>0, (3.37) 

we conclude (A-9F(0)) _1 G £(/i 1+5 (S x ), /i 4+ *(S 1 )) for Re A > A* as in the proof of Lemma EH 
The second step consists of proving the resolvent estimate 

, -1 



|A| (X-dF(0)Y 



< c , Re A > A» . 



Since Kcq n < for \n\ sufficiently large, elementary calculations show that 

|A-<?„| 2 > c |A| 2 , Re A > A» , neZ, 
for some Co > 0. Thus, setting := A(A — q„) _1 it follows 

sup \S*\ < 00 . 

nGZ\{0} ,Re A>A, 

Similarly, there is c\ > such that 

1-^ — Qn\ 2 > c\\q n \ 2 , RcA>A*, neZ, 

and we thus obtain from 



(3.38) 

(3.39) 
(3.40) 

(3.41) 



n+l I 



A - q n 



Qn+l — qn 



I)' 
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together with (j3~B"4"|) . ^jD| , and ([3"Hj) combined with (|53T|) that 

sup |n| \S* +1 - S*\ < oo . 

nGZ\{0},ReA>A» 

Finally, noticing that the right hand side of 
Qn+2 — 2<?„+i + q„ 



(3.42) 



< \z. 



n+2 1 



H(n + 2) - 2p(n + 1) + fi(n) 



/i(n + l) 



- Zn+l| « 2 + 



|z„ - z n+ i\n 2 



is bounded by (|3.32p and writing 
71,2 \Sn+2 ~ 2-5„_|_i + <S , ^ +2 | 



| | \ n3 Qn+l 



A - g„ 



|A(g n+2 - 2q n+ i + q n ) + q n (q n +i - q n +2) + q n +2(q n +i - q n )\ 



we deduce from ([331]) . (|3~35]) . (|3~39l) . and (pUP]) 



sup n 2 1 5^ +2 - 2S* +1 + S* | < oo 



(3.43) 



nGZ\{0} , Re A>A* 

Therefore, Lemma EH] and (j3T40|) . (|Q2|) . and (|3T43|) imply 



|A| (X-dF(0)) 1 



< c , Re A > A* 



□ 



£(C 1 +' 5 (S 1 ) 

whence (|3.38|) due to (|3.37[) . This proves the assertion 

Now we are in a position to establish a wcll-posedncss result regarding equation (|3.28[) . 

Theorem 3.11. There exists an open zero neighborhood O C V in h i+s (§ 1 ) such that for each 
Po G O there is T :~ T(po) > an d a unique solution 

p E C([0,T},V) DC 1 ^^},^^^ 1 )) 

to 

d t p = F(p), t>0, p(0) = po- 

Moreover, p([0, T])) C O. 

Proof. We shall invoke [H Thm.8.1]. Fix £ 6 (0, S) and put ■& := (5 - 0/3. Set 



(S 1 



in [14j Thm.8.1]. As 8 <E (0, 1) was arbitrary in Proposition 13. 101 it follows that 

-dF(0) G ^(/^(S 1 ),/^^ 1 )) . 

Thus, since H(/i 4+ '»(S 1 ), /i 1+ ^(S 1 )) is open in C(h i+ ^ (S 1 ), /i 1+ ^(§ 1 )) , there is an open zero neigh- 
borhood C/ c in /i 4 +«(S x ) such that -<9F(p) G H(/i 4+ «(S 1 ), (S 1 )) for each p G U f . Then 
:= (~1 W with W from Lemma T3.8I is an open zero neighborhood in /i 4+,5 (§ 1 ). Furthermore, 
dF(p) : /^(S 1 ) /^(S 1 ) for p G O is the part of -dF(p) G W (/i 4+ « (S 1 ) , ^(S 1 )) in 

^(S^-^+^S^^^+^S 1 ))^ 

with continuous interpolation functor (•, ^ and 

{/i G ^(S 1 ) ; dF(p)[h] G /^(S 1 )} = /^(S 1 ) . 

Now the assertion is a consequence of [H] Thm.8.1]. □ 
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To finish off the proof of Theorem 12.11 let 

pGC([0,T],V)nC 1 ([0,T],/ l 1+5 (§ 1 )) 

be the solution to (|3.28[) for a given initial value po G O. Then 

Qi := S{p,d tP ) e C([0,T] ) h a + a (S i )) , Qo ~T{p,Q- l -K,{ P )) e C([0,T]> 2+5 (S°)) 

by Lemma f3.3[ Proposition 13.41 and Proposition 13.61 Since p solves (|3.12|) . it follows, for p = 
with t € [0, T] fixed, that 



ld t p -1 



a du„(ct>PQ ) + p dr p (<f> P *Qo)) on T 



Recalling from ([33]) that A(0£Q o ) = and div (zx V(</>£Q )) = in we deduce from Gauss' 
Theorem 

77^ / (a Q d Vp (tfQo) + Mr p {€Qo))da = — L [ (a d v (<%Q ) + Mr (€Qo) ) da 
l w o| Jr p l w o| 

and thus, due to 4>*Q = Qo and (|3.1ip on RS 1 , that 



|ViVp| 



der = 0, te[0,T]. (3.44) 



Consequently, with 

W) : = 0£ (t) QiW e ^ + ^ (t) ) , p (t) := 0S w Qo(*)/i 2+ *(n° w ) 

for £ <G [0, T] we obtain a solution (p, P;, P a ) to (|2.11[) - (|2.14[) which is unique up to additive constants 
in the pressures Pi and P a . This yields Theorem 12. II 

4. Proof of Theorem 12.21 

We first prove instability of the trivial solution if g- t > g a as claimed in Theorem 12.21 Recall 
that jj := Qjio 2 /2. 

Theorem 4.1. If Qi > g D , then 

d t p = F(p) , t > , p(0) = pa , 

has backward solutions which do exponentially decay to zero. In particular, the trivial solution 
p = of this flow is unstable. 

Proof. The compact embedding /i 4+l5 (§ 1 ) ^ /i 1+<5 (§ 1 ) and Proposition l3~10l imply that dF(0) has 
compact resolvent. So the spectrum of dF(0) consists of eigenvalues only, which, according to 
Proposition 13.91 are given by {q n ; n £ Z \ {0}} and q := 0. Since Qi > g , (|3.29j) implies that 
p(l) > 0, hence (|3.30j) shows that Regi > 0. Clearly, 

inf {Re q n ; Re q n > 0} > . 

The assertion now follows from [HI Thm.9.1.3]. □ 

To prove stability of the trivial solution if p Q > Qi, we need an auxiliary result. 
Lemma 4.2. Let W be given as in Lemma \S.8[ Then 

[ (1 + p)F{p) da = , peW. 
Js 1 
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Proof. Fix p G W and set 



fcoiWH/e^CS 1 ); / (l + p)/d<r = 



' G 1 

We claim that 1 — IZ(p) is an isomorphism on /^^(S 1 ). To see this, set Tf := T(p, S(p, /)) for 
/ e /^(S 1 ). Then, as in ([37H|) . 

/ (l + p)ft(p)/da = / (l + / 9)B o (p)T / da = 

and Lemma 13.81 implies that 1 — IZ(p) is indeed an isomorphism on /ij"^ 5 (S 1 ) . But, as above, 

B {p)T(p,-JC{p)) ehl+'iS 1 ) . 

Therefore, 

= (1 -^(p))- 1 B„(p)T(p J -JC(p)) e ftJ+^S 1 ) . 

□ 

We conclude the proof of Theorem 1 2. 2 1 by stating the stability result for which we need to dchnc 

^(S 1 ) := {/e/^S 1 ); f /d<r = 0}, 

for s > 0. 

Theorem 4.3. 7/ £> > Pi, the trivial solution p = of 

d t p = F{p) , t > , 

is stable. More precisely, there are numbers oj, r, M > suc/i that for each initial datum po G /ij +(5 (§ 1 ) 
with \\po\\h 1 + s (s 1 ) ^ r fiene is a unique global solution 

p G C(R+, ^o +5 ( §1 )) n G 1 (S 1 )) 

wit/i p(0) = po and 

P(t)(l + ^f] +||p(t)(l + p(t))||, ll+ , (sl) < Me- 4 p„(l + ^) t>0. 

Proof. Letting £ := p + p 2 /2 for /? G W, problem c\ p = F(p), t > 0, is equivalent to 

a t C = G(C), t >0, (4.1) 
where G G C 2 (Z,h 1+S (S 1 )) with Z := {C =p + p 2 /2; p G W} is given by 



G(C) := (2V1 + C- l)F(2v^+C-l) • 
Moreover, Lemma BJ implies G € G 2 (Z , /iJ+^S 1 )) for Z := Z n h^ +s (S 1 ). Also, 

dG(0) = dF(0) G ^(^(S 1 ),^^ 1 )) • 

Thus (9G(0) has compact resolvent and its (point) spectrum equals {q n ; n G Z \ {0}}. Now, by 
(1^291) . 

while Q > Q\ implies that 

p{n) = \n\ {a - 2( 7o - 7i ) - an 2 ) < -2( 7o - 7i ) < , n G Z \ {0} . 

This, combined with (|3.30[) . shows that the spectrum of <9G(0) is contained in a half plane 
[Re A < — lj] for some ui > 0. The assertion now follows from [HI Thru. 9. 1.2] applied to (14. ip . □ 
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Note that our analysis yields an explicit estimate of the exponential decay rate uj in terms of 
the physical parameters through (|3.29p . 
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